Convergence of the two-point function of the 
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Abstract We consider the two-point function of the totally asymmetric simple ex- 
clusion process with stationary initial conditions. The two-point function can be 
expressed as the discrete Laplacian of the variance of the associated height func- 
tion. The limit of the distribution function of the appropriately scaled height func- 
tion was obtained previously by Ferrari and Spohn. In this paper we show that the 
convergence can be improved to the convergence of moments. This implies the con- 
vergence of the two-point function in a weak sense along the near-characteristic 
direction as time tends to infinity, thereby confirming the conjecture in the paper of 
Ferrari and Spohn. 



1 Introduction and result 



The totally asymmetric simple exclusion process (TASEP) is arguably the simplest 
non-reversible interacting stochastic particle system, and it is also one of the most 
studied. Particles live on Z and they satisfy the exclusion constraint: each site can be 
occupied by at most one particle. Therefore a particle configuration can be denoted 
by tj € {0, 1} Z , where rj ; - = means that site j is empty while tj ; - = 1 means that 
the site is occupied. The dynamics of the TASEP is then defined as follows: every 

Jinho Baik 

Department of Mathematics, University of Michigan, Ann Arbor, MI, 48109, USA. 
e-mail: baik@umich.edu 

Patrik L. Ferrari 

Institute for Applied Mathematics, Bonn University, Endenicher Allee 60, 531 15 Bonn, Germany, 
e-mail: ferrari@uni-bonn.de 

Sandrine Peche 

U.F.R. de Mathematiques, Universite Paris Diderot, 75205 Paris, France, 
e-mail: peche @ math . univ-paris-diderot . fr 



1 



2 



Jinho Baik, Patrik L. Ferrari, and Sandrine Peche 



particle tries to jump to its right neighbor with rate one. The jumps occurs only if 
the exclusion constraint is satisfied. 

It is known [12] that the only translation-invariant stationary measures of the 
TASEP are Bernoulli product measures with parameter p £ [0, 1], namely, 

F(77/ = 1) = p forally'eZ. (1) 

Here p is the average density of particles. The cases p = and p = 1 are trivial and 
in the following we fix p e (0, 1). This system is referred as stationary TASEP. 
The two-point function is defined as 

5(;,r):=E(7 7 . / (0'?o(0))-p 2 . (2) 

Note that this equals the covariance of T]j{t) and 77o(0). Hence the two-point func- 
tion carries the information on how site j at time t is correlated with site at time 0. 
It is known that 

Y t S(j,t)=p(l-p)=:X (3) 

jez 

and also S(j,t) > 0. This implies that jS(j,t) can be thought of as a probability 
mass function in j £ Z. Indeed this equals the probability that a second class particle, 
which was at site at time 0, is at site j at time t [9]. It is also known that the 
expectation of j with respect to the probability mass function j^S(j,t) satisfies 

E^ = (1-2P>, (4) 

jez * 

and the variance scales as [16, 18] 

£_.2^£)_ ((1 _ 2pW 2 = 0( ,4/3 ) _ (5) 

jez X 

as t — > oo. Therefore, for large time t , one expects the scaling form for S as 1 



16 6sc V 2j 1 / 3 f 2 / 3 ) 2j 1 /3f2/3 

for some non-random function g sc . The precise expression of g sc was first conjec- 
tured in [15] based on the work [6]: 

g sc (w) = [ s 2 dF w (s) (7) 
Jr 

where F w (s) is the distribution function defined (17) below. 

In order to understand the presence of the second derivative in (6) and the second 
moment formula (7), we recall that TASEP can also be seen as a stochastic growth 



1 The multiplicative factor A was incorrectly written as \ in [14]. This is a typographical error. 
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interface model, whose discrete gradient of the height equals 1 — 277. The dynamical 
rule is that when a particle jumps to the right, a valley \/ changes to a mountain 
/\. More precisely, let N t (j) denote the number of particles which have jumped 
from site j to j+ 1 during the time interval [0,f], and define the height function 



ht(j) 



2M(0)+I/ =1 (1-2tj ( -(0) forj>l, 

2N t (0) for 7=0, (8) 

2N t (0) -Z? =y+1 (l -277,(0) for; < -1. 



Then initially ho(0) — and ho(j) — ho(j — 1) = 1 — 2t7j(0), and at the instance a 
particle jumps from site j to j+l, the height function at position j increases by two. 
Note that h t (j) — ho(j) = 2N t (j). It was shown in [14] that the two-point function 
can be expressed as 

S0\0 = t(4Var(M-)))C/) (9) 

with A being the discrete Laplacian, (A f)(j) — f(j - 1) —2f(j) + f(j+ 1). Since 
it is known that F w (s) has mean [6], this explains the presence of the second 
derivative in the conjectured formula (6) and the second moment formula (7). 

Define the probability distribution functions of the location-rescaled height func- 
tion, 

F w (s,t) :=P(M[(1 -2p)t + 2w X l/3 t 2 / 3 }) 

> (1 - 2%)t + 2w(l - 2p)j 1 / 3 f 2 / 3 - sx 2/3 t : / 3 ) . (10) 

The function F w in (7) (which is defined in (17) below) was conjectured in [15] to 
be the limit 

]imF w {s,t)=F w {s). (11) 

The convergence (11) for each s was later proved in [10]. This strongly indicates 
the validity of (6). A missing part in concluding (6) is the convergence of the mo- 
ments of F w (s,t) which is a stronger statement than (11). Our main result is that the 
moments indeed converge. 

Theorem 1. For all £eN, 

lim / s e dF w (s,t) = [ s e dF w (s) (12) 
uniformly for w in a compact subset ofR. 

As a consequence we obtain the convergence of the two-point function is a weak 
sense. 

Corollary 2 We have, with % := p (1 — p), 

lim2 Z 1 /V/ 3 5([(l -2p)t + 2wx il3 t 2 l\t) = ^' s '» (13) 
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if integrated over smooth functions in w with compact support. 

The proof of this corollary is given in Section 5. An improvement of the analysis 
in this paper can yield the convergence in the point- wise sense in (13). However, we 
do not consider this direction in this paper. 

For completeness, let us state a formula of the limiting distribution function F w (s) 
explicitly. Let P u be the orthogonal projector on the interval [u, +°°). Set 

KMs(x,y) ■= / Ai(x + s + X)Ai(y + s + X)dX, 

JR+ (14) 

FgueM := det(l -P Q K A ^P Q ). 

Fgue is the GUE Tracy- Widom distribution function [17]. We also define the func- 
tion 

g(s,w):=e-^( f e w{x+ ^Ai(x+y + s)dxdy + f % s (x)p s (x,y)$ w , s (y)dxdy ) , 

(15) 

where 

$ w ,s(x)--= f e wz+ws K M , s {z,x)dz, %. s (x) ■■= I e wz Ai(x + z + s)dz, (16) 

JR- JR- 

and p s (x,y) := (1 -PoK AUs P )- l (x,y). Now 

F w (s) := |- (F GVE (s + w 2 )g(s + w 2 ,w)) . (17) 

There is an alternative formula expressed in terms the Lax pair equations of the 
Painleve II equation obtained in [6]. But we will only use the formula (17) in this 
paper. One can also consider the joint distributions for different values of w and a 
formula can be found in [5]. 
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2 Setting and strategy of the proof 



The height function h t (j) associated to a TASEP with any initial condition can be 
related to the last passage time of a directed last passage percolation (DLPP) model. 
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Over the last decade or so, the so-called solvable, or determinantal DLPP models 
[8, 1 1, 13] were studied extensively. These are the models for which the probability 
distribution of the last passage time can be expressed explicitly in terms of Fredholm 
determinants. The DLPP model corresponding to the stationary TASEP is not one 
of solvable models but can be related to one after suitable analytic continuation of 
the parameters. This yields the following formula of F w (s,t). 
Fix w e R. Let us set 2 (recall that J = p (1 - p )) 

2m = (l-2x)t + 2w(l-2p)x l/3 t 2 ^, 2d = (1 - 2p)t + 2wx l/3 t 2 ? 3 , (18) 
and define the functions 3 

- e a(x-y) f ( +0 \m-d 

L{x,y) = — — A e -^-y) } Z ^ P > dz for x > y, 

K 2m. J T (l-p-z) m+d 

p (19) 

e a{x-y) r (i_ _ \m+d 

R( x ,y) = — — J> e z ^ \ P . ' . dz for x < y, 
y ,yi 2m Jr_ p (p+z) m - d yi 

with 

a:=\-p. (20) 

We define the kernel 

K m4 {x,y) = [ L{x,z)R{z,y)dz, (21) 

and the distribution function 

F(u) :=det(l- P u K m4 P u ). (22) 



Finally, we set 
where 



Go(«) =gi(u)+g 2 (u)+g 3 (u), (23) 



2ad — m 



«!(»)=»+ 1/4 _ fl2 . 
g2(u) = (\lfa,PuK m ,dY-a), ^ 

g 3 (u) = (K; n /t-P u )y a A^-PuK m4 P u )- l P u {l-K m4 )y- a ), 



2 To be precise, we need to take the integer parts of the formulas since m and d need to be integers. 
Since the error between the formula above and the integer parts is 0(1), this does not result in 
any significant changes in the estimates and hence for convenience we define m and d as in (18) 
without restricting them to be integers in this paper. However, we remark that if we restrict m and 
d to be integers, one occasionally need to be careful in the precise formulation of the estimates and 
the exposition becomes more involved. We do not discuss these subtleties in this paper. 

3 For any set of points 5, the notation § Fs f(z) dz denotes the integral over a simple closed contour 
which encloses the points S but excludes any other poles of the function /. The contour is oriented 
counter-clockwise. 
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with y a {x) = e~ ax - Then it was shown in [10] that 4 

M*,*) = ( f/ y)i/3 ^ (F(u(s,t))Go(u(s,t))) (25) 



where 
Set 



u = u(s,t) :=t + sX~ 1/3 t 1/3 . (26) 



G (s,t):=G (u(s,t)), F(s,t):=F(u(s,t)). (27) 

The main technical part of this paper is the following estimates: 5 

Proposition 1 (Uniform upper tail estimates). There exist positive constants so, to, 
c and C such that 



s-F(s,t)G {s,t)/(t/x) 



1/3 



<Ce~ c ^, s>s Q , t>t Q . (28) 



The bound holds uniformly for w in a compact subset o/R. 

Proposition 2 (Uniform lower tail estimates). There exist so, to, c and C such that 

\F w {s,t)\<Ce- c ^ V \ s<-s , t>t Q . (29) 

The bound holds uniformly for w in a compact subset o/R. 

Theorem 1 now follows. 

Proof of Theorem 1. We only consider I > 2. The case £ = 1 follows easily. We first 
write the integral on the left-hand-side of (12) as the sum of the integral over R + 
and the integral over R_. For the integral over R + , integrating by parts twice and 
using the fact that F w (-,t) is a cumulative distribution function, 

/ /dF w (s,t) =-£(£- 1) / J-2( s -F( S ,t)^&h)ds (30) 

JVL+ M + \ {t/X) /J 

for £ > 2. It was in [10] that in addition to (11) we also have limit 
F(s,t) Fg\je(s + w 2 )g(s + w 2 1 w) for each s as t — > °°. Thus due to Propo- 

sition 1 the Lebesgue dominated convergence theorem can be applied and we find 
that (30) converges to 

-£{£-\)J^ s ( - 2 (s-FavE(s + w 2 ) 8 (s + w 2 ,w)ys. (31) 



4 The formula (25) is the formula (4.10) of [10] when b = -a if we take into account (26) . See 
(5.21) of [10] for the formula of the function G (u) = G a - a (u). 

5 The exponents of the bounds are not optimal. The bound in (28) and (29) can be improved to 
£ g -c|s| 3 / 2 Ce~ c l s l 3 respectively. The improved bound for (28) can be achieved if we keep track 
of a slightly better estimate in the analysis presented in this paper. On the other hand, in order to 
improve the bound (29), we need a different approach such as Riemann-Hilbert analysis as in [3,4]. 
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On the other hand, integrating by parts once, 

/ s e dF w (s,t) =-£ s e - l F w (s,t)ds. (32) 

JR- JR- 

Thus again, using the Lebesgue dominated convergence theorem can be applied due 
to Proposition 2 and from (1 1) we find that (32) converges to 

-t [ s £ - l F w {s)ds. (33) 

JR- 

Integrating (31) and (33) by parts backwards and using the fact that F w is a cu- 
mulative distribution function, we find that the sum of these two integrals is the 
right-hand-side of (12). 

□ 

The estimate (28) for the upper tail is obtained by analyzing the formulas (22) 
and (23) asymptotically using the saddle-point analysis. This asymptotic analysis is 
very close to that of many previous papers, for example [2, 10, 1 1]. We use some of 
the results directly or improve upon them. See Section 3. 

For the estimate (29) on the lower tail, we note the following. Consider the 
TASEP with step-initial condition i.e. rj ; -(0) = 1 for j < and 77, (0) = for j > 1. 
Then the associated height function hf ep (j) satisfies h,Q £p (j) = \j\. This means that 
initially ho is bounded above by /iQ tep . Since the initial condition of the stationary 
TASEP is independent of the dynamics, we find that h t is stochastically bounded 
above 6 by hf ep . Hence 7 

P(^0')>«)<^f ep (;')>«)• (34) 

But P(/zf ep (j) > u) is known to be precisely F{u) of (22) [11]. Therefore we have 
F w (s,t) < F(s,t) = det(l -P u K m4 P u ). (35) 

Thus the estimate (29) follows if we show that F(s,t ) is bounded above by Ce~ c > s > 
for negative large enough s. This in turn follows if we show the same bound for the 
Fredholm determinant (22). For this purpose we follow the idea of Widom [19] 
which seems not as well-known as it should be. See Section 4. 

6 This can also be seen easily from the corresponding directed last passage percolation (DLPP) 
models. The DLPP model for the stationary TASEP is the DLPP model for the TASEP with the 
step initial condition plus an extra row and an extra column with non-zero weights. 

7 We would like to thank Ivan Corwin and Eric Cator for communicating this observation with us. 
This observation simplified the proof of the lower tail estimate which we originally obtained by 
estimating F w (s,t) directly. 
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3 Proof of Proposition 1: upper tail 



The proposition follows from (39), (37) and (51) in the below. 



3.1 Asymptotics for F 

The function F(u) = det(l — P u K m ^P H ) is the distribution function of the last pas- 
sage time of the directed last passage model with i.i.d. exponential random vari- 
ables. It is well-known [11] that this also equals the distribution function of the 
largest eigenvalue of the Laguerre unitary ensemble (LUE) which is defined as 
M m ,d = j^XX* where X is a (m — d) x (m + d) random matrix with i.i.d. standard 
complex Gaussian entries. This equality can also be seen explicitly in Appendix C 
of [10] where K m( j was shown to be same as the correlation kernel of the LUE up to 
a conjugation by a multiplication. The asymptotics of LUE and F(s,t) = F(u(s,t)) 
were considered in several papers, especially in [2, 10, 1 1]. We have: 

Lemma 1. Fix so e R. Then 

limF{s,t)=F GVE {s + w 2 ) (36) 

t— 

uniformly for s £ [so, 00 ) an d w in a compact subset ofR. Furthermore, for given 
io £ K and to > 0, there exist positive constants C and c such that 

\l-F(s,t)\ <Ce~ cs (37) 

for s > sq and t > to- 

The bound (37) can be found in, for example, Section 3.1 of [2]. 8 



3.2 Evaluation of g\ 

A direct computation using (18), (20), and (26) shows that 9 

*i(«) = »+ 1/4^2 ' • (38) 

8 The exponent of the upper bound is not optimal: the optimal exponent is e~ c l s l 3/2 . But we do not 
consider such an issue in this paper. 

9 The formula becomes s(t/x) 1 ^ + 0(1) where 0(1) is independent of s if we take the integer 
parts in the definition of m and d in (18). This is an example of the subtleties mentioned in the 
Footnote 2. This results in the additional term 0(t ~ 1,/3 ) in (39). Since this is not a function in s, we 
cannot obtain the bound (CI). However, this issue can be fixed by shifting s to s — 0(1) /(t/%) 1 ^- 
In other words, the centering and scaling u = t + s(t j%) '/ 3 needs to be changed slightly to reflect 
the difference of the formula of (18) and their integer counter-parts. 
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This implies that 

- i|r = -^"•') £ w^ +*<> - *M> <»» 

The term 1 — F(s,t) can be estimated using (37) and F(s,t) is bounded by 1 since 
it is a distribution function. We now show that gi(u) / (t / %) x l 3 and gi(u) / {t / xY ^ 
are uniformly (in t ) bounded by exponentially decaying functions in s. 

In the rest of this section, we only consider the case when w > 0. If w < 0, 
we need to start with a different decomposition of Gq(u) ((5.22) instead of (5.21) 
of [10]). After this change, the analysis is completely analogous. For the case when 
w = 0, we can proceed as in the case when w > but with a yet slight modification: 
see (6.31)-(6.34) of [10]. We skip the detail when w < and w = 0, and assume 
from now on that w > 0. 



3.3 Estimations on g2 and gs 



Recall the definition (24) of g2(u). It is a direct calculation to show that (see (3.15) 
of [10]) 



R{x,y)y-a{y)dy = Z(p)y- a {x), Z(p) :-- 



(l-p) m+d 

r>m—d ' 



(40) 



forxe M, for a e (-1/2, 1/2). Using this, g2{u) = {y a ,P u L{i -Po)V-a)- Inserting 
the formula \j/ a and L(x,y), we obtain 



S2{u) = / <s%(x+y)dxdy, 



where 



*):= 



2ni 



-z(u+x) _ 



(z + p) 



m—d 



jdz- 



(41) 



(42) 



(1 -p-z) m+d 
Thus (see (6.19) of [10]) 

(t/xr 1/3 g2(u) - f H t {x+y)dxdy, H,(y) := (f/z) 1/3 ^Wf/l) 1/3 ). (43) 

Similarly, recall the definition (24) of ^(m). Using (40), an argument similar to 
that for (41) implies that 



(l-K m . d )yf- a (x)=e a 



1 - / ,3%(-u+x + y)dy 



(44) 



We also note that, similar to (40), we have (see (3.15) of [10]) 
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Wa{y)L(y,x)dy 
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Wa{x) 



(45) 



forx e R, for a e (-1/2,1 /2). Using this, we find that 

K; nJ (i-p u ) Wa ( x ) 

[ M{-u+x+y)dy- [ J%(-u+x+y)J%(z + y)dzdy 



= e 



where 



J%(x) :-- 



1 



e^ i 7 p -* m : d d , 



2%\Z{p)J T _ p ~ (z + p) m - d 

This implies that we can express (see (6.26)-(6.28) in [10]) 

(f/xr 1/3 s3(«) = {^,At% 

where 



(46) 



(47) 



(48) 



H t (y + ^ )dy - H t (x + y)H t (y + ^ )dxdy 

1- / H t (y + $)dy 

m + 



(49) 



with H t (y) = (t/x) 1/3 ^(y(t/x) 1/3 ) and H t (y) = (f/z) 1/3 ^(y(f /l) 1 / 3 ), and the 
operator A t is defined by A t = Pn(l — K t )~ 1 Po where the kernel of K t is 



K t {^ 2 )=e w ^~^ j H t (x + ^)H t (x+^ 2 )dx, &>&>0, 

m + 



(50) 



and K t , £2) = otherwise. 

We obtain the following estimates for g2 and #3. 

Lemma 2. 77iere are positive constants c and C such that 



(t/xr 1/3 82(u) 



<Ce~ 



\(t/x)- lli 83{u)\<Ce- 



(51) 



for all s > and t > 0. 

Proof of Lemma 2. Note from the formula (42) that J^(x) = J%(x;u) is a function 
ofx + M . Hence // r (y) —H t (y;s) is a function of y + 5. Thus, //, (y;s) =H t (y + s;0). 
The same holds for H t (y) = H t (y; s). 

Basic bounds for the functions H t (y) and H t (y) were obtained in (6.15) of [10]: 
for any j3 > there exist positive constants Cp and such that 



\H t {y;s)\<C p e-Py and \H t (y;s)\ < C 1 ^ 



(52) 
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uniformly for t > 0, y > 0, and s > 0. In particular, the bound holds for H, (y; 0) 
and H t (y ; 0) when s = 0, for t > and y > 0. Thus using H t (>■; s) = (y + s; 0) and 
inserting y + s in place of y in (52), we find that: for any j3 > there are positive 
constants Cp and C'^ such that 

\H t (y;s)\<C p e-^ and |ff ( (y;j)| < C^"^ (53) 

uniformly in t > 0, y > and 5 > 0. 

The bound for (f/z) _1/3 ^2(«) follows from (43) and (53). 

We now estimate \(t / %y l ^g3{u)\- Choosing j3 > \w\, (53) implies that 
\<P t (%)\ <Ce-$ s e-W- w ^ for a positive constant C. Thus, 

\m L 2 { M +) <C'e-P s , (54) 

for a constant C' uniformly in t > and s > 0. On the other hand, (53) implies that 
\%(^)e w ^ | is bounded by a constant. Since we assume w > (see Section 3.2), we 
find that H^H^p^) is uniformly bounded in t > and s > 0. Finally, using the 
inequality 

I \A t || < 1 1(1 - 1 1 + 11(1- ^w+J- 1 - (1 - K,)- 1 1| (55) 

where K Ai w 2 +s is the Airy kernel restricted on (w 2 + s,°°), and the fact (see (6.36) 
of [ 1 0] ) that 1 1 ( 1 - K Ai ) - 1 - ( 1 - K, ) - 1 1 1 -)• as t ->• °o imply that 1 1 A t \ \ is uniformly 
bounded in t > and 5 > 0. Therefore, the bound for (t/%) ~ l ^g?,(u) follows from 
\(<P t ,A t ¥ t )\ < ||0 f ||||A,||||^||. 

□ 



4 Proof of Proposition 2: lower tail 

Recall from Section 3.1 that K m ^ is a similarity transform of the correlation kernel 
of the LUE M m d . Since the correlation kernel of the LUE is a positive projection, 
all the eigenvalues, which we denote by )J,j,j = 0, 1,2,- • • , of P u K md P u are real and 
/I, e [0, 1]. It was shown in Appendix B.3 of [10] that [ij e [0, 1) if u > 0. From 
this we find that det(l - P u K m . d P u ) = Uj>o(l - My) < Uj>oe-^ = e - Tr ^ K ^" p -\ 
Therefore, 

F(s,t) < exp(-Tv (P u K m . d P u )). (56) 

This trick is due to Widom [19]. 

The trace has the following lower bound: 

Proposition 3. There exist positive constants to, sq, c such that 



Tr(P u K m . d P u )>c\s\V 2 



(57) 
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for all s < —sq and t > to- 

The same estimate was obtained in the context of random permutations and an 
oriented digital boiling model by Widom [19]. We follow the paper [19] to prove 
the Proposition, and as such we only sketch the main ideas and do not provide all 
the details of the proof. Once this proposition is proved, then Proposition 2 follows 
from (35) and (56). 

Proof of Proposition 3. Since the operator^, is trace class with continuous kernel, 
we have 

\m+d -tn—d 



r -1 r r e wu (1 



(l_ w )m+rf ft-d dwdz 



z) m+d w m - d (w-z) 2 

(5o) 

-1 r f e MF «^ dwdz 



(27ri) 2 7r, Jr e MF ^ (w-z) 2 

where 

F„(z):=ii , z-lnz + yln(l-z), u' := ^. (59) 

Here M:=m-d and y:= ^ = (^) 2 + C^f- 1 / 3 ). Note that u' = jj + O^ 1 / 3 ) 
if s is in a bounded set, and u' > for all s < 0. We analyze (58) asymptotically 
using the saddle-point analysis. Note the presence of the singularity ^j^p in the 
integrand. 

We first consider the case where 

(i-Vr) 2 + e<« < (i + vr) 2 -sof~ 2 / 3 , (60) 

for some e > (small, but fixed) and sq » 1 also fixed. The critical points are F u 
are 

^ {u ' ] = u ~ ± ir L± h^ {u ' (1 + ^ )2)(M ' = (1 = ^ )2) - (61) 

The two critical points are non-real and |z*(«')l = < p < 1. Consider the fol- 
lowing two contours: 

wHz+|e ie , O<0<27z;, (62) 

and 

Z= l + |z+-l|e ie , O<0<2&, (63) 

respectively. Then 



k + | 2 kl 2 



(64) 
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Thus along these contours Re (F u ) achieves its relative maximum (resp. minimum) 
at zf. Hence these paths are of steep-ascent and steep-descent for F u . We chose 
to work with these explicit contours instead of the contours of steepest-ascent and 
steepest-descent for convenience. Due to this reason, we need to modify the contours 
locally near the critical points if u' is close to (1 + ^Jj) 2 . Namely, in this case, the 
contours above become almost tangential and are almost parallel to the direction 
along which Re(F u ) is constant. Then we cannot apply the saddle-point method. In 
this case, we simply modify the contours locally near the critical points so that it 
passes through the critical points along the steepest descent direction as pictured in 
Figure 1 for the z-contour. A similar modification is needed for the w-contour. This 
small modification does not yield any significant changes in the estimation. For the 
convenience of presentation, we work with the above explicit contours and skip the 
details on how the formulas changes after the modifications. The same procedure 
was also explained in Section 6.2 of [7] for the similar estimations. 

We now deform the original contours in (58) to the new contours of steepest- 
ascent and steepest-descent, which we call by the same names, Jo and 1]. We first 
deform the original contours to those in (a) of Figure 2 where fo is the contour of 
steepest-ascent and the part of n except for the segment from z~ to z+ is the part of 
the contour of steepest-descent. These contours can be divided as in (b) of Figure 2 
and we have 

-1 f f e MF »M dzdw 1 r f e Mf »W dzdw 

( ) ~ (2^i]2 RV - fr /r, (w-z) 2 + J^tif 2 h %(' (w-z) 2 ' 

(65) 

Here the first integral needs to be interpreted as the Principal Value due to the di- 
vergent terms in the integrand. The second integral is from the contributions of the 
pole in the deformation of the contours. The contours in the second double integral 
are defined as follows. The w-contour, ^ ', is a segment from z^ to z[ to the left of 1 
and to the right of 0. The z-contour, I]", encircles the whole segment ^ but not 1, 
see Figure 2. 
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Fig. 2 The subdivision of the integration from (a) the ones in (58) to (b) the ones in (65). 



Setting Q(z) := exp(MF H (z)), the Cauchy's integral formula implies that the sec- 
ond integral of (65) equals 

-1 fQ»_ dw = -M(F u{z i)-F u{z -))_ (66) 



2m }%' Q(w) 2m 



Noting that F u (z+) = F u (z c ), we have 
.m) 2 !<gfr(' 



e MU^MU. dzdw ^ Mlm[Fu{zt) y 



(2m) 2 JvTr;' (w-z) 

Observe that when u 1 = (1 + y/f) 2 , the two critical points coincide and 
Z c := z^r = T^fi- In addition, F M{ y + ^i(z c ) & M. Thus 

-Im(F H (z+)) = lm(F M{l+V7)2 {z c )) -Im(F„(z+)) 

= [ ^-FM{v))dv. 

Using the definition (59) of F u and the fact that z+(v) is a critical value, we find that 
^F v (zt(v)) = }}zt{v). From the formula (61) of z+ 

Im(z+(v'M)) = ^"^"^ ((1 + Vff v') 1/2 

zv (69) 

>|((i + Vr) 2 -V) 1/2 

since v' satisfies the condition (60). Therefore, (68) implies that 

(67)>M^((l+Vr) 2 -«') 3/2 - (70) 

Recall that ^7 = (1 - p)/p + O (f~ 1/3 ), M = p 2 t(l + 0(f~ 1/3 )), and u' = u/M with 
u = t + s(t/x) 1 ^- Then, we can choose a sq > large enough (but fixed indepen- 
dently of t) such that for all s < —sq it holds (1 + y/j) 2 — u'> — c\st~ 2 l 3 for some 
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c\ > 0. Therefore for u' satisfying (60), there is a positive constant c such that 

1 r r eMF„(w)-MF u (z) 

7^2 f —( v—dzdw>c{-sfl 2 (71) 

(2n\) 1 J'g Jr{' (w-zY 

uniformly in t . 

We now show that the contribution of the Principal Value integral in (65) is much 
smaller than (71). Indeed we will show that this is 0(1). This proves (57) by taking 
the constant c smaller than one in (71). 

A direct computation shows that 

F"(z + ) = (1 - M-TTyr^c-^ (?2) 
u{Zc) { J) (zt) 2 (z+-l) 2 ■ { ' 

This implies 

\F:(zt)\~((l+VY) 2 -u') 1/2 ~s 1/2 t- l/3 (73) 
as u' ->• (1 + yy) 2 - st ~ 2 / 3 , while for (1 - ^/f) 2 + e < u' < (1 + ^/y) 2 - e we have 
\K'(zt)\ = °( 1 )- Thus > for ^ general u' satisfying (60), cif" 1 / 3 < |F H "(z+)| < c 2 
for some positive constants c\ and c^. Hence 0(t 1 ^) < \/MF"(zt) < 0(t l l 2 ). Let 
us choose the parts Vo(zf) and Vi(zf) of the paths fo and I] respectively whose 
size are B(MF"(z^))~ 1 / 2 . Then both parts become smaller as t — > °°. Because 
the paths fo and T\ are chosen to be steep descent, the contribution coming from 
roxri\{Vo(z+)UV (z-)}x{Vi(z+)UVi(z 6 7)} is at most of order 0(1) if B is 
chosen large enough (but fixed). Let us first consider the contributions from the 
intersecting contours Vb(z^) x Vi(z+) and Vo(z^) x Vi(z~ ). Due to the symmetry, it 
is enough to consider the contribution of Vo(z^ ) x V\ (z+), given by 

ll 7 +\ • = _1 PV [ [ e MF u {w)-MF u (z) \ dz dw (74) 

Now we have to see if this integral is bounded by a constant. Since z converges to 
z+, we use the Taylor's series of F u in z — z£. Since z^ is a critical point, the function 
F u (z) in the exponent may be approximated as F„(z+) + 5^'(zc)(z — ) 2 - It can be 
checked that the contributions from the higher order terms are negligible. Changing 
the variables as z = z+ +z'(M J Fj'(z+))- 1 / 2 , w = z+ + w'(M J F i ;'(z+))- 1 / 2 , we obtain 

(2ra) 2 JiKJM (w'-z') 2 

which is finite. 

Let us now show that the contribution of the non-intersecting contours 
Vo(zf) x Vi(zJ) are also bounded from above by some constant. To that aim, B 
being fixed, we assume that so is chosen large enough so that sq 3> B. This time the 
singularity term l/|w — z| 2 is bounded from above and one can easily deduce that 

for all u' e ((1 - Vr) 2 + e, (1 + Vr) 2 - ^ 2/3 ), 
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Lp.y. / / e MF u ( W )-MF u (z) _}_ dzdw < 0(1) . 

I 2 Jv n (z±) M(zJ) (w-z) 2 



Combining the whole, we have shown that the contributions from the first integral 
in (65) is 0(1) and (57) is proved for u ! e [(1 - V?) 2 + £, (1 + V?) 2 - s t ~ 2 / 3 ). 
We now consider the case where 

u'e (o,(i-vr) 2 + £)- (76) 

In (60), we could have chosen e > small enough so that 

(i-Vr) 2 + £<^((i + Vr) 2 + (i-Vr) 2 )- (77) 

Consider the Laguerre Unitary Ensemble j^XX* where X is a (m — <i) x (m + 
random matrix with i.i.d. complex standard Gaussian entries. Denote by 
> A2 > • • • > Kn-d its ordered eigenvalues. By the definition of the correlation 
kernels, we have 

(m— d \ 
I 1/(A,)J , (78) 

where / = («', +°°). This can be bounded below as 

(m—d \ /m-d \ 

E 1/(A,)J > E (I l/ e Wj, (79) 

where I £ = ((1 — ^/y) 2 + e,+°°). Now, we call on the results of [1], giving con- 
vergence rates for the spectral distribution of random sample covariance matrices. 
Let F m _^ denote the empirical probability distribution function associated to the 
spectral measure: 

i m— d 

F,n-d(x) = I U,< x - (80) 

m a j=l 

Let also F be the cumulative distribution function of the Marchenko-Pastur distri- 
bution p defined by the density 

dp J (u t , —x)(x — u c ) , 

JL = vv + A -±\ w u c, (x , (81) 

dx 2kx " 

where y = and uf = (1 ± ^/y) 2 . It is well known that F m _ rf (x) F(x) a.s. for 
all x. In [1] it is proven that 

max|E(F m _ d (jc))-F(x)| < (m-d)-^ 2 . (82) 

jc>0 

Then (78) and (79) imply that 

Tr(K t )> {m- d){\ - - jyf + e)) - {m - d) 1 ' 2 . (83) 
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With the condition (77) on e, F((l — y 7 /) 2 + e) < 1 uniformly and since 
m — d = p 2 t + 0(t 2 / 3 ) -> oo, we find that there exists a positive constant C = C(e) 
such that 

Tr{K,)>C{m-d) (84) 

uniformly in t and for u' satisfying (76). Now as u = t — s(t/%) 1 / 3 and u > 0, we 
have {-sfl 2 = (t - ufl 2 (% jtfl 2 < % l l 2 t. Thus since m — d = p 2 t + 0{t 2 l 3 ), (84) 
implies that there exists a positive constant c such that 

Tr (K t ) >c{-s) 3/2 (85) 

uniformly in t and for u' satisfying (76). Thus (57) is proved for u' satisfying (76). 
This completes the proof of Proposition 3. 

□ 



5 Proof of Corollary 2 



Let us consider the rescaled height function 



„. , h t {j{w)) - [(1 -2 X )t + 2w(l -2p) X W] 

H,(W).= _ z 2/3 f l/3 ' (86) 

with j(w) = (1 - 2p)f + 2wx 1 / 3 f 2 / 3 . By (10), F w (s,t) = F(H,(w) < s). We have: 



(87) 



G,(w) :=War(H,(w)) = [ s 2 dF w {s,t) - ( f sdF w (s,t)) , 



and, in the original variables, 

Var(h,U(w))) = (x 2/3 t^) 2 G t (w). (88) 
Using the notation 8 := (2j 1 / 3 f 2 / 3 )- 1 , by (9) 

/ M 2,VV/3^ M| , )/(wMw = X ^A»^)-2 G ,W + G,(w-^) /M<i>| , 

16 jk o z 

(89) 

By Theorem 1 and the fact that f K 5c/f w (i) = (see [6]), we have that G t (w) con- 
verges to g S c(w) uniformly for w in a compact set of E. Therefore, for smooth test 
functions / with compact support, as t — > °° this expression converges to 

4 / 5sc(w)/"(w)Jw= ^ / ^'»/(w)Jw. (90) 
16 Jk 16 Jr 
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